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1. 


INTRODUCTION 


In  flight  dynamic  studies  of  aircraft,  there  is  frequently 
a  requirement  for  the  estimation  of  aerodynamic  derivatives.  Theoretical 
approaches  to  this  estimation  often  require  a  conformal  mapping  from  a 
unit  circle  to  the  aerofoil  shape  concerned  as  one  of  the  steps  involved. 
Therefore  an  investigation  has  been  carried  out  into  the  transformation 
from  circles  into  practical  shapes.  Some  results  of  this  investigation 
are  presented  in  this  Technical  Memorandum,  which  is  essentially  a 
compendium  of  shapes  (aerodynamic  and  non-aerodynamic)  in  the  z-plane 
obtainable  by  conformal  mappings  from  the  unit  circle  in  a  £ -plane. 

Shapes  obtained  by  varying  the  position  and  size  of  this  circle  are 
too  numerous  to  be  catalogued  and  are  considered  to  be  essentially  just 
variants  of  the  main  ones  considered  here. 

The  method  of  conformal  mapping  from  the  exterior  of  the  unit 
circle  to  the  exterior  of  a  desired  closed  curve  has  widespread 
applications,  for  example  In  aerodynamics  and  in  solid  mechanics. 

However,  comparatively  few  shapes  have  been  specifically  given  in 
applications.  The  best  known  are  the  ellipse  and  Joukowski  aerofoils, 
which  are  variants  of  the  narrow  slit,  and  the  Karman-Tref f tz  aerofoils 
which  are  variants  of  areas  bounded  by  two  circular  arcs  with  equal 
curvature  (see  Kober  [1]  and  Glauert  [2])  and  hypocyloids  (see  the  recent 
paper  by  Clement  [3]).  Most  other  shapes  are  generated  only  by  multistep 
and  ad-hoc  mappings  (see  Clement  (4]  and  Wu  [5],  [6])  or  by  infinite 
series  (see  Abbott  and  Doenhoff’s  book  [7]  on  the  method  by  Theodorsen 
and  the  technique  by  Naiman) .  The  present  work  aims  to  bring  to  light 
another  range  of  curves,  besides  the  well-known  ones,  which  can  be 
produced  by  simple  analytic  formulae  or  by  simple  numerical  methods. 

These  curves  are: 

1.  Double-cusped  aerofoils, 

2.  Biconvex-aerofoils  with  different  leading  and 
trailing  edge  angles, 

3.  Biconvex-aerofoils,  by  integration, 

4.  Triangles  with  arbitrary  vertex  angles. 
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5.  Regular  polygons, 

6 .  Rhombs , 

7.  Rectangles, 

8.  Aerofoils  by  finite  product  formulae. 


The  first  two  types  of  aerofoil  curve  are  derived  from 
synthesised  formulae.  The  next  five  types  of  curve  are  obtained  from 
the  Integration  of  dt/d(.  The  last  type  of  aerofoil  curve  is  obtained 
from  a  finite  product  formula  subjected  to  some  constraint  as  developed 
previously  in  Ref.  9. 


The  following  sections  will  treat  those  types  of  curves  and 
their  generating  methods  in  the  above  order.  Of  particular  concern  is 
the  closure  condition  for  the  integration  method.  This  is  automatically 
satisfied  for  Schwarz-Christof f el  transformations  but  not  necessarily 
for  the  type  of  transformations  used  to  generate  the  shapes  from  3.  to 
7.  A  discussion  of  this  condition  has  not  been  found  in  literature 
available  to  the  author.  Section  6  describes  the  new  integral  method 
used  to  plot  the  curves  in  the  transformed  region  and  the  Fourier 
integral  method  used  to  evaluate  the  complex  coefficients  of  the  descending 
series,  which  are  given  in  the  corresponding  figure,  for  each  transformation. 
A  novel  way  to  generate  an  approximate  square  is  also  mentiqned  there. 


By  the  use  of  these  transformations,  two-dimensional  problems 
of  pure  circulation,  inviscid,  incompressible  flow  about  these  shapes 
can  be  easily  solved.  Similarly,  they  can  be  used  to  solve  the  problem 
of  incompressible  Stokes'  flow  about  these  shapes.  The  two-dimensional 
problem  of  stress-concentration  for  an  infinite  plate  can  also  be  solved 
for  these  shapes  by  the  use  of  the  well  known  complex-variable  method. 


2.  DOUBLE-CUSPED  AEROFOILS 

The  formula  for  this  aerofoil  was  synthesised  by  closely 
following  Joukowskl's  formula.  Wo  guess  a formtil n  of  the  fora 


iSNjdx  j  in  tv, ,  ifl  >:.)  iv  aaonao-i.mi 
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and  select  a,  b  and  k  such  that  dz/d(  vanishes  only  at  £  -  ±  1  and  to 
zeros  of  order  one.  The  resulting  formula  Is 


*(c)  *  C  + 


(l-a)a  C 
(1+a)  a  ’ 


which  has  Its  derivative  as 


dz  (1-a 

d?  1+a 


(1-a)2  ;2+a 

1+a  (C2-a)! 


This  Is  clearly  a  function  with  the  desired  properties.  The 

constant  a  determines  the  maximum  thickness  of  the  aerofoil  Figure  1 

gives  such  an  aerofoil  with  the  constant  a  being  0.07.  The  image  of  a 

nest  of  concentric  circles  is  plotted  using  a  graphic  plotter.  The 

small  irregularities  on  these  curves  are  due  to  the  limitation  on  the 

resolution  of  the  plotter  as  well  as  of  the  integration  subroutine 

used.  The  small  table  in  the  same  figure  gives  the  values  of  the 

coefficients  a  for  the  descending  series 
n 


*(c)  -  a, 5  +  E  *  C 

l  _  i  n 

n"-l 

which  is  a  Laurent  expansion  of  the  transformation  function  z(c). 
Transformations  for  aerofoils  in  subsequent  figures  also  have  their 
series  given  in  this  way. 

The  double-cusped  aerofoils  with  different  leading  and 
trailing  halves  are  given  by  the  formula 


1  .  (1-a) ‘ 


Its  derivative  is 


dz  _  ,  .  2b_  1  (1-a)2  ;2+2b;-t-a 

dc  ‘  +  1+a  ?  ‘  (1+a)  *  (C2-a)2 


(4) 


which  vanishes  only  at  £  *  Si  to  zeros  of  order  one.  The  constant 
a  determines  the  maximum  thickness  of  the  aerofoil  and  the  constant 
b  determines  the  extent  of  asymmetry  between  the  leading  and  trailing 
edge.  Figure  2  is  a  plot  of  an  aerofoil  so  generated  with  a  ■  0.07, 
b  -  -0.3. 


3.  AEROFOILS  WITH  FINITE  AND  DIFFERENT  LEADING  AND  TRAILING 
EDGE  ANGLES  (GENERALISED  KARMAN-TREFFTZ  AEROFOILS) 


From  the  Karman-Trefftz  formula 


z  "  u 


,  0<pS2, 


we  construct  the  new  formula 
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(5) 


for  which  the  derivative  vanishes  only  at  t  ■  ±1  to  zeros  of  order 
(u— 1 )  and  (Ap-l)  respectively. 


In  the  above  formula,  u  is  the  external  angle  for  the  trailing 
edge,  Au  is  the  external  angle  for  the  leading  edge  (0SA£l),  k 
affects  the  ratio  between  the  influence  by  the  leading  edge  and  that  by 
the  trailing  edge  (0<<<a‘)  ,  and  n  affects  the  propagation  of  the 
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circular  arcs  from  the  leading  and  trailing  edges  (0<n<l).  For  the 
special  case  where  X  is  equal  to  one,  the  above  formula  reduces  to 
the  standard  Karman-Tref f tz  formula. 

Figure  3  gives  such  an  aerofoil  with  u  “  1.9,  A  *  0.7, 
ic  «*  0.4  and  n  -  0.4.  Effects  of  each  individual  parameter  on  the  shape 
of  the  aerofoil  are  demonstrated  in  figures  4,  5  and  6.  In  these 
latter  figures,  individual  parameters  ic,  n,  X,  are  set  to  1.0,  0.2, 

0.6  respectively  while  the  remaining  three  parameters  are  left  the 
same  as  in  figure  3. 

Analysis  easily  shows  that  the  curve  given  by  equation  (5) 
asymptotes  circular  arcs  of  the  Karman-Tref ftz  type  at  the  leading 
and  trailing  edges  of  the  aerofoil.  By  using  all  four  parameters 
at  our  disposal,  the  basic  NACA  0015  aerofoil  with  no  camber  is 
approximated  with  the  choice  p  ■  1.9,  X  ■  0.5263,  <  m  0.22  and  n  *  0.53 
The  resulting  aerofoil  is  given  in  figure  7.  Note  that  this  aerofoil 
is  the  image  of  the  unit  circle  in  the  C -plane.  The  images  of  other 
circles  in  the  c  plane  can  give  even  better  approximation  of  the  NACA 
aerofoil. 


4.  SHAPES  DERIVED  FROM  THE  INTEGRATION  OF  dz/di;  IN  THE 
FORM  OF  FINITE  PRODUCTS 


It  has  been  known  (see  Carrier,  Krook  and  Pearson  [8] 
that  the  integration  of 


dz 


i  n  2y 

tt  n  (;-$.)  t  | ♦.  |  *  i» 

4  k-i 


n 

E 

k=l 


1 


gives  a  transformation  from  the  outside  of  a  unit  circle  to  the 
outside  of  a  polygon  with  some  choice  of  {4*^}.  Here  we  point  out 
that  transformation  by  integration  of  dz/d c  need  not  be  a  haphazard 
process.  Indeed,  if  we  have 
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and 


n 

Z 

k-1 


Vk 


0 


(7) 


then  the  integration  of  dz/d £  will  give  a  single-valued  function  z(c) 

for  |^|  >  1,  such  that  its  derivative  is  given  by  the  first  member 

of  equation  (6).  The  closure  equation  (7)  has  not  been  found  in 

available  literature.  This  condition  here  is  arrived  at  by  requiring 

that  z(0  is  single-valued  and  by  considering  the  contour  integral  of 

the  first  member  of  equation  (6)  along  a  circle  of  very  large  radius. 

The  vanishing  of  this  integral  has  (7)  as  it  necessary  and  sufficient 

condition.  (Another  result  from  the  condition  of  closure  is  that  if 

Y  “  <Y, *  Yt . .  )  i*  a  vector  of  real  values  of  external  angles  of 

.  1  i  n 

a  polygonal  then  equation 

n  ix 

T  y  e  «  0  , 

k-1 

will  have  at  least  a  real  vector  solution  x  -  (x.  ,  x_,  ...,  x  ).  Each 

.  l  L  n 

shape  of  the  polygon  has  a  unique  solution  x  -  (x^,  x^.  . . . >  xn) »  all 
Xj  within  [0,  2ir],  determined  up  to  a  constant  common  to  all  components. 
Proof  is  by  the  Rienann  mapping  theorem  for  the  exterior  of  the  polygon 
from  the  region  |;|  >  1). 

Our  application  of  the  transformation  (6)  subject  to  the 
condition  (7)  gives  the  biconvex  aerofoil  of  figure  8.  The  derivative 
dz/d z,  is  taken  as 
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with  a  ■  0.8,  8  ■  0.6,  c  ■  0.142.  The  condition  (7)  is  satisfied 
with  this  choice  of  parameters.  Note  that  a  determines  the  leading 
edge  angle,  6  the  trailing  edge  angle  and  c  is  a  function  of  a  and  g 
by  equation  (7).  This  way  of  generating  aerofoils  gives  more  flexibility 
than  the  use  of  formula  (3). 

The  generation  of  a  triangle  with  given  vertex  angles  is  also 
a  special  case  of  equation  (6)  subject  to  the  closure  condition  (7). 

The  derivative  for  this  case  is 


dz 

d; 


,  10.2-a-g.  ia.a.  ib.g 

(4~e  )  (c-e  )  (;-e  ) 


Figure  9  is  a  triangle  corresponding  to  a  •  0.5,  3  =  0.7,  a  »  2.094 
and  b  -  3.807.  Note  that  a  and  b  are  chosen  such  that  equation  (7) 
is  satisfied  with  the  values  of  a  and  8  given. 

/ 

Shapes  with  some  kind  of  symmetry  usually  have  condition 
(7)  automatically  satisfied,  and  are  given  in  text  books  (such  as  [8]) 
without  any  mention  of  the  latter  condition.  The  following  curves 
belong  to  this  category. 


By  choosing 


dz 


1 

4 


(en-i) 


2/n 


we  have  an  n-sided  polygon.  The  examples  are  given  in  figures  10, 

11,  12  and  13.  As  n  tends  to  infinity,  the  polygon  becomes  a  circle. 
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Aerofoil  shapes  can  also  be  generated  by  the  use  of  finite 
product  formulae  as  considered  in  a  previous  study  [9].  This  method 
consists  of  using  formulae  such  as 

r(0~zc  -  (c-l)1’9(C+0.1+0.2i)0'1(C-0.1+0.2i)'1  , 

with  the  constant  zc  chosen  such  that  z-C  is  of  order  C  1  or  smaller 
as  |c|  tends  to  infinity.  A  computer  program  is  used  to  ensure  that 
all  singularities  of  dz/d{  are  contained  in  the  unit  disc  jcjsl  with 
one  at  t;«l  to  generate  a  finite  acute  angle  at  the  trailing  edge. 

(The  computer  programs  for  checking  the  singularities  of  dz/dt  and 
the  initial  shapes  of  these  aerofoils  were  written  by  Mr.  C.  A.  Martin 
The  author  acknowledges  his  help  to  sustantiate  this  method).  An 
aerofoil  generated  by  this  method  is  given  in  figure  16.  The  constant 
i.  has  been  calculated  in  the  computer  program  and  the  value  of  z(0 
is  plotted  in  the  figure. 
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6.  PLOTTING  OF  CURVES  AND  COMPUTATION  OF  DESCENDING  SERIES 

To  plot  curves  in  the  transformed  region  we  use  the  formula 
developed  in  Ref.  [9] 

z(C)  “  /  el2lTtda(t)  +  exp [/  log(5-el2rt)  [da(t)+dt]  ] 

0  0 

where  a(t)+t  is  given  by 

27ri[a(t)+t]  «  log[z  (ei2lTt)  -  z^],  ObtSl  , 

i2rt 

and  Zp  is  an  arbitrarily  chosen  point  interior  to  the  closed  curve  z(e  ). 
This  method  of  plotting  was  used  to  generate  the  curves  in  all  figures 
of  this  Memo.  The  advantage  of  this  method  is  that  infinitesimal 
features  close  to  the  unit  circle  are  preserved,  in  contrast  to  the 
descending  power  method  where  a  very  large  number  of  terms  are  required 
close  to  the  points  of  singularity  of  dz/d£  and  this  number  becomes 
infinite  as  any  singularity  is  approached. 

Far  away  from  the  aerofoil  it  is  quite  economical  to  use 
decending  series,  as  demonstrated  by  the  shapes  of  the  images  of 
concentric  circles  in  Figures  1  to  17.  The  descending  series  for  each 
aerofoil  is  given  in  the  table  in  its  corresponding  figure  as  described 
previously  in  section  2.  The  complex  coefficients  of  these  series  are 
computed  using  Fourier  integrals  on  z(e^2ltt).  This  method  is  applicable 
here  as  the  image  z(c)  is  known  for  each  point  r  =  e*2lTt  Gf  the  unit 
circle  in  the  ^ -plane. 

As  a  point  of  curiosity,  we  note  that  the  descending  series 
for  the  function 

z  “  x  +  iy 
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Two-dimensional  problems  of  incompressible  Stokes'  flow  about 

these  cross-|’3ections  and  of  sti ^ss+concentration  for  an  infinite 

<  <' 

elastic  plate  with  a  hole  of  these  shapes  can  also  be  solved  using 
the  conformal  mappings  given  here.  £  c, L  • ,  , 

t-7  ri  o5;  -t>f  0  •  'C 


/ 


M  U*  J  ANjmr  M)*j  XV  U  J  Jl  ILK  n  1. 1  m 


10  - 


where 


< 


< 


V. 


cos(2nt)  if  jcos(2irt)|  S  — 

7i 


COs(2TTt) _ 

I  cos(2irt)  |  Jl 


if  I  cos(2itt)  |  > 


sin(2nt)  if  1 sin(2nt)|  £  — 

v2 


sin(2fft) 


I  sin(27rt)  I  /2 


if  I  sin(2-nt)  I  > 


/2 


1_ 


gives  an  approximate  square,  as  in  figure  17.  The  reason  is  that  the 
coefficients  of  its  ascending  series  are  mostly  zero  or  very  small 
and  decrease  in  magnitude  very  rapidly  with  increasing  power  of  z. 

It  is  also  convenient  to  note  in  this  Memo  that  the  n-corner  starred 
shape  used  in  some  stress  concentration  problems  is  generated  by  the 
formula 

1/n 

z  m  (^n  +  +  1  +  m)  ,  / 

cn 

where  n  is  positive  integer  and  m  is  a  positive  constant  less  than 
unity. 


7.  CONCLUSIONS 

.  .  /,  *  r  •  ;•/  /U  1  I  .  1  <•  •  ' 

It  has  been  shown  here  that  there  are  many  practical^shapes 
obtainable  from  conformal  mapping  using  simple  techniques.  Employing 
the  results  given  here,  it  is  only. a  matter  of  st might  forward 
application  to  calculate  the  pure  circulation,  inviscid,  incompressible 
flow  about  the  aerofoils  of  this  Memo.  The  pressure  distribution  then 
follows  directly. 
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